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Abstract. The results of the quantitative investigations of the renormal- 
ization of the absorption edge of different compounds by the isotope effect are 
described. The obtained dependence of the band gap energy versus isotope 
mass is in reasonable agreement with estimated value obtained from the tem- 
perature dependence of E ff in LiH crystals. For the first time it was shown that 
the dependence E g = f(M) for the different materials has a parabolic charac- 
ter ln{dE g /dM) = 6.105(lnE g ) 2 - 7.870(lnE g ) + 0.565. Nonlinear character of 
this dependence may be indicated on the nonlinear dependence of the potential 
scattering on the isotope mass. 

PACS: 32.10.B + 71.20.P + 78.20 + 78.40.F 

Isotope substitution is a well-defined and easily controllable method to inves- 
tigate intrinsic renormalization mechanisms of elementary excitations of solids. 
First of all, phonon frequencies are directly affected by changes of the average 
mass of the whole crystal or its sublattices (in the case of virtual crystal ap- 
proximation) even when we look upon them as noninteracting particles, i.e. as 
harmonic oscillators. The renormalization of the fundamental electronic gap 
by electron-phonon interaction also depends on the isotope mass. Measuring 
the energy gaps in samples with different isotopic composition then yields the 
difference in the changes of the valence- and conduction-band renormalization. 

In this communication we report the first results of the quantitative study the 
dependence of the band-gap energy on the isotope effect for different compounds. 
As was mentioned above isotopic substitution only affects the wavefunctions of 
phonons; therefore, the energy values of electron levels in the Schrodinger equa- 
tion ought to have remained the same. This, however, is not so, since isotopic 
substitution modifies not only the phonon spectrum, but also the constant of 
electron-phonon interaction. It is for this reason the energy values of purely 
electron transitions in molecules of hydride and deuteride are found to be dif- 
ferent (see, e.g. [1]). This effect is even more prominent when we are dealing 
with a solid [2]. This conclusion was confirmed on a qualitative level as early 
as the 1930s in Ref. [3]. Intercomparison of absorption spectra for thin films of 
LiH and LiD at room temperature revealed [3] that the longwave maximum (as 
we now know to be the exciton peak (see, e.g. [2])) moves 64.5 meV towards 
the shorter wavelengths when H is replaced by D. As will be shown below, this 
effect becomes even more pronounced at low temperatures. 
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The mirror reflection spectra of mixed and pure LiD crystals cleaved in liquid 
helium are presented in figure 1. For comparison, on the same diagram we have 
also plotted the reflection spectra of LiH crystals with clean surfaces. All spectra 
have been measured with the same apparatus under the same conditions. As 
the deuterium concentration increases, the long-wave maximum broadens and 
shifts towards shorter wavelengths. As can clearly be seen in figure 1, all spectra 
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exhibit a similar long-wave structure. This allows us to attribute this structure 
to the excitation of the ground (Is) and the first excited (2s) exciton states. 
The energy values of exciton maxima for pure and mixed crystals at 2K are 
presented in Table I. The binding energies of excitons calculated by the 
hydrogen-like formula, and the energies of interband transitions E 9 (details see 
[4]) are also given in Table I. 

The ionization energy, found from the temperature quenching of the peak of 
reflection spectrum [4] of the 2s state in LiD is 12 mcV. This value agrees fairly 
well with the value of AE 2s calculated by the hydrogen-like formula. Moreover, 
Eb = 52 meV for LiD agrees well with the energy of activation for thermal 
quenching of frec-exciton luminescence in these crystals [2] . 

At the weak scattering potential the mean-square vibrational amplitude 
<u 2 > of an atom dependes on the phonon frequencies and the eigenvectors, 
the atomic masses (at low temperatures), as well as the temperature (at high 
temperature). Isotope substitution results in a slightly different vibrational 
amplitudes (epecially at low temperature) and phonon frequencies and in first 
approximation equivalent to changing the temperature (see also [5]). The mass 
dependence of <u 2 > becomes vanishingly small at temperature on the order 
or higher than the Debye temperature. Changes in either the isotope masses 
or temperature thus lead to changes in the band gap (see above Fig. 1) via 
the electron-phonon interaction, even at zero temperature in the case of mass 
changes. In general, the renormalization of the band gap [4] and its temperature 
dependence result from a complicated interplay of 1) first- and second-order 
electron-phonon interactions that contribute to the energy of conduction and 
valence bands; 2) changes due to thermal or isotopic lattice expansion, and 3) 
changes in the phonon occupation numbers (see [5] and references therein). 

The dependence of the band gap energy on isotopic composition has already 
been observed for insulators and lowest (indirect-direct) gap of different semi- 
conductors (see also [6]). It has been shown to result primarily from the effect 
of the average isotopic mass on the electron-phonon interaction, with a smaller 
contribution from the change in lattice constant. This simplest approximation, 
in which crystals of mixed isotopic composition are treated as crystals of iden- 
tical atoms having the average isotopic mass is reffcred to as virtual crystal ap- 
proximation (VCA) [7]. Going beyond the VCA, in isotopically mixed crystals 
one would also expect local fluctuations in the band gap energy from statistical 
fluctuations in local isotopic composition within some effective volume, such as 
that of an exciton. As follows from Fig. 1 excitons in LiHxDi-^ crystals display 
a unimodal character, which facilitates the interpretation of their concentration 
dependence. Fig. 2 shows the concentration dependence of the energy of inter- 
band transition E ff . As can be seen from Fig. 2 VCA method cannot describe 
observed experimental results. As will shown below this deviation from linear 
low (VCA approximation) is connected with isotope-induced-disorder in isotope 
mixed crystals LiH^Di-a;. 

The temperature and isotopic mas dependence of a given energy gap E 9 (T,Mj) 
can be described by average Bose-Einsteinm statistical factor n# corresponding 
to an average phonon frequency 9i as (see also [5; 8] 
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E 3 (T,M,) = E fcar - a r (^f )V2 [i + 2 n B ], (1) 

where = l/[exp(^) -1] and E;, ar and a r the unrenormalized (bare) gap 
and the renormalization parameter, respectively. In the low-temperature limit, 
T«9i, equation (1) reduces 

E s (T,Mj) = E fcar - a r (^) 1/2 (2) 

Here E fl (T,Mj) is independent of temperature and proportional to (1/Mi) 1//2 , 
whereas a r «s the energy difference between the unrenormalized gap (Mi — ► oo) 
and the renormalized value [5] . 

In the high-temperature limit, T>>8i and Eq. (1) can be written as 

E s (T,Mj) = E bar - 2T^, (3) 

and E s (T,Mj) is independent of Mj [5]. The extrapolation of Eq. (3) to 
T = OK can be used to determine the unrenormalized gap energy E{, ar , i.e., 
the value that corresponds to atoms in fixed lattice position without vibrations 
(frozen lattice [4]), from the measured temperature dependence of E ff (T) in the 
high-temperature (i.e. linear in T) region. 

Using Eq. (2) it can be written the difference in energy AE 9 between a given 
energy gap in isotopically pure material (LiH) and its isotope analogue (LiD) 



AE 9 = E 9 (M 4 ) - E g (M nat ) 



Mi 



(4) 



As can be seen from Table 1 and results of [5] AE g at 2K equals AE g = 
0.103 eV and E g (LiH, T = OK) =5.004 eV (linear approximation and E g (LiH, 
T = 300K) = 4.905 eV then using Eq. (3) we obtain a r = 0.196 eV. This 
magnitude is very close (approximately 84%) to the value of 0.235 eV of zero 
vibration renormalization of the energy band gap in LiH crystals (details see 
[4;5]). Using Eq. (4) we obtain AE g (theor) = 0.134 eV that is very close, on 
the other hand, to observed experimental value equals 0.103 eV. The discrepancy 
between these values may be caused by the negligible contribution of the isotopic 
lattice expansion to the band gap renormalization. 

The single-mode nature of the exciton reflection spectra of mixed crystals 
LiHrDi — x agrees qualitatively with the results obtained with the virtual crystal 
model (see, e.g. [7; 9]), being at the same time its extreme realization, since 
the difference between ionization potentials (£) for this compounds is zero. Ac- 
cording to the virtual crystal model, ( = implies that AE g = 0, which in 
contradiction with the experimental results for LiH^Di-a; crystals (details see, 
also [6]). In the light of obtained dependence E g = f(T, Mj) for LiH^Di^a, it is 
very interested to study such dependence for different materials. Although the 
details of the renormalization process of E 9 are not known, observation of the 
isotope renormalization of E g for many materials force to pay a great attention 
on this effect. 

By now the change in E g caused by isotopic substitution has been observed 
for many broad-gap and narrow-gap compounds. Below we briefly discuss the 
variation of the electronic gap (E 9 ) of different crystals with its isotopic com- 
position. In the last decade the whole row of different semiconducting crystals 
was grown. These crystals are diamond [10], copper halides [11; 12], germanium 
[13;14], Si [8] and GaAs [15]. All enumerated crystals show the dependence of 
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the electronic gap on the isotope masses. It should be noted that the indicated 
effect (the variation of E s and Eh, see Table I and II) have maximum values in 
LiH crystals, although this effect in other crystals with isotopic composition are 
currently being reliably measured and investigated well. All numerated crystals 
show the dependence of the electronic gap on the isotope masses. 

Before we complete the analysis of these results we should mention that 
before these investigations, studies were carried out on the isotopic effect on 
exciton states for a whole range of crystals by Kreingol'd and coworkers [16- 
21]. First of all we should name the classic crystals CU2O [17-19] with the 
substitution O 16 — * 18 and Cu 63 — > Cu 65 . Moreover, there have been some 
detailed investigations of the isotopic effect on ZnO crystals [6] where E g was 
seen to increase by 55 cm" 1 (O 16 -> O 18 ) and 12 cm" 1 ( at Zn 64 -> Zn 68 ) 
[16;20]. In the paper of Krcingol'd and coworkers [21] it was shown that the 
substitution of a heavy S 34 isotope for a light S 32 isotope in CdS crystals resulted 
in a decrease in the exciton Rydberg constant (E(, ). 

More detailed investigations of the exciton reflectance spectrum in CdS crys- 
tals were done by Zhang et al [22] . Zhang et al to directly determine the binding 
energy (E&) and the corresponding band gaps (E g ) from a hydrogenlike model. 

In this manner it was obtain binding energies of A excitons of 26.4±0.02 meV 
and 26.8±0.02 meV in 112 CdS and ™ a *CdS, respectively. The corresponding 
band gaps E 3 (r 7c - T 9v ) arc 2.5806(2) cV and 2.5809(2) eV, respectively. In 
the case of B excitons, these values are Ef, = 27.1±0.2 meV (27.1±0.2 meV) 
and E 3 (r 7c - T 7v ) = 2.5964(2) [2.5963(2)cV] for the 112 CdS [™ a *CdS] sample. 
Unfortuntealy, the n =2 excited states of the A and B excitons could not be 
observed in other isotopic CdS samples. Better samples are required for such 
measurements . 

For GaAs or ZnSe, isotope substituents of either type should lead to shifts of 
the Eo gap which have been calculated to be 430 (420) and 310 (300) /ieV/amu 
for cation (anion) mass replacement, respectively [15;23]. These values are in 
reasonable agreement with data measured for GaAs [dE /dM Ga = 390 (60) fieV/amu] 
and preliminary results for isotopic ZnSe obtained by Zhang et al [22] based on 
photoluminescence measurements of the bound exciton (neutral acceptor Ii) 
[dE/dMse = 140 ± 40 /ieV/amu and dE/dM Zn = 240 ± 40 /j,eV/&mu] . 

Such behavior, however, is not found in wurtzite CdS. A previous reflectivity 
and photoluminescence study of ™ a *Cd 32 S and " at Cd 34 S shows [21] that for 
anion isotope substitution the ground state (n = 1) energies of both A and B 
excitons have a positive energy shift with the large rate of dE/dMs = 740±100 
/xeV / amu. This value is more than one order of magnitude larger than <9E / dMcd 
obtained by Zhang et al [22]. 

Several groups have conducted low-temperature studies of the direct and 
indirect band gaps of natural and isotopically controlled Ge single crystals. For 
the first time Agekyan et al [13] used photoluminescence, infrared absorption, 
and Raman spectroscopy with a Ge crystals of natural composition and a crys- 
tals with 85 % 76 Ge and 15% 74 Ge. They found an indirect band-gap change 
AE S = 0.9 meV and a direct band-gap change AE 9 = 1.25 meV with an error 
of ±0.05 meV. Etchegoin et al [24] and Davies et al [25] reported photolumi- 
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nescence studies of natural and several highly enriched , high quality single 
crystals of Ge. Measurement of the energies of impurity-bound excitons per- 
mits by Davies et al the direct determination of band-gap shifts with the crystal 
isotope mass because the radiative recombination does not require phonon par- 
ticipation. As may be expected from the very large Bohr orbit of the excitons 
(see Davies et al [25] and reference therein), their binding energy only depends 
on the average isotope mass and not on the isotopic disorder (see, however [4]). 
The rate of band-gap energy change with isotope mass as determined by Davies 
et al is dE /G /dA = dE NP /dA = 0.35±0.02 mcV/amu. Etchegoin et al [24] 
obtained a very similar value. 

Zollner et al.[25] have performed a numerical calculation of the electronic 
bands using an empirical pseudopotential method including the necessary lat- 
tice dynamics. They found for Ge (dEja/dA) e ^ p = 0.41 meV. The total cal- 
culated shift of the indirect band-gap energy with isotope mass adds up to 
(dE/G/dA)i 0ta ; = 0.48 meV. This result compares favorably with the experi- 
mental values stated above by Davies et al and by Etchegoin et al who reported 
(dE /G /dA) tota , = 0.37±0.01 meV/amu. 

Measurements of the direct band gap at the T point (k — 0) in the Brillouin 
zone have also been performed. Though the direct band gap is technologically 
less important than the minimum indirect band gap, determing the dependence 
of this gap on isotope mass is of the same fundamental significance as the indi- 
rect band-gap studies. Davies et al [27] used low-temperature optical-absorption 
measurements of very thin samples of Ge single crystals with natural composi- 
tion and three different, highly enriched isotopes. They found dE/dA = 0.49 
± 0.03 meV/amu for the temperature extrapolated to zero. Parks et al [28] 
have used piezo- and photomodulated reflectivity spectra of four monoisotopic 
and one natural Ge crystals. These techniques do not require the extreme sam- 
ple thinning which is necessary for optical-absorption measurements and the 
derivative nature of the spectra emphasizes the small changes. The excellent 
signal-to- noise ratio and the superb spectral resolution allowed a very accu- 
rate determination of the dependence of Edg on isotopic mass. At very low 
temperatures an inverse square-root dependence accurately describes band-gap 
dependence E DG = Eg G + ^= . 

A fit through five data points yields E^ G = 959 meV and C = -606 meV/amu 1 / 2 . 
Written as a linear dependence for the small range of isotopic masses, Parks et 
al find dE DG /dA = 0.49 meV/amu, in perfect agreement with the results of 
Davies et al [27]. Parks et al also determined the isotope mass dependence of 
the sum of the direct gap and the split-off valence band (A ) and find d(E£> G 
+ A )/dA = 0.74 meV/amu. The experimental results can be compared to 
the Zollner et al [26] calculations which arc found to be of the correct order of 
magnitude. The theoretical estimates for the contributions of the linear isotope 
shifts of the minimum, indirect gaps which are caused by electron-phonon in- 
teraction, are too large by a factor of ~ 1.7 and for the smallest direct gap they 
arc too large by a factor ^3.2. 

Substitution of Ga 70 on the Ga 76 increases the band gap in GaAs [23] on 



5 



10.5 cm -1 . The interested results were communicated in papers of Cardona and 
coworkers [15; 12], where it was studied the dependence of E g on the isotope 
effect in CuCl crystals. When the Cu 64 on the Cu 65 is substituted the value 
of E 9 in CuCl crystals decreased on 1.24 cm -1 , e.g. the isotope effect on the 
electronic excitation has an opposite sign. 

Considering the series of Ge, GaAs, ZnSe, CuBr, for example, the 3d states of 
the first constituent play an increasing role in determining the band structure. 
In Ge these states can be considered as localized core states (atomic energy 
level w -30 eV). Already, however, in GaAs they have moved up in energy by 
10 eV, and their hybridization with the top of the valence band affects the gap 
. Proceeding further in the series, this effect becomes more important, and 
in CuBr and Cu 3d states even overlap in energy with halogen p-states, with 
which they strongly hybridize. Therefore, we cannot exclude that the main 
reason for the opposite sign of the isotopic effect in these compounds may be 
connected to the different character of the d-electron-phonon interactions in 
these semiconductors [6]. 

Analogous investigations of the direct absorption edge of Si (E g = 3.562 eV) 
have been performed in paper [8] made the value of the coefficient 9E ff /i9M = 2.0 
meV/amu. The measurements of dE g / dM for isotope effect in CsH (CsH^CsD) 
gave the value <9E ff /<9M = 60 meV/amu [29]. The change of the indirect gap of 
diamond between pure C 12 and C 13 has been determined by Collins ct al [10], 
using for this purpose the cathodoluminescence spectra of diamond. From the 
results of Collins et al it was concluded that the dominant contribution arises 
from electron-phonon coupling, and that there is a smaller contribution due to 
a change in volume of the unit cell produced by changing the isotope. These 
two terms were calculated as 13.5 ± 2.0 and 3.0 ± 1.3 meV respectively (details 
see, also [30]). 

All of these results are documented in Table II, where the variation of E ff and 
dEg/dM are shown at the isotope effect. We should highlighted here that the 
most prominent isotope effect is observed in LiH crystals, where the dependence 
of E;, = f (Ch) is also observed and investigated (see, also [2]). Using the least- 
squares method it was found the empirical dependence of ln<9E 3 /<9M ~ f(lnE 9 ), 
which is depicted on Fig. 3. As can be seen the indicated dependence has a 
parabolic character: 

ln(dE g /dM) = 6.105(lnE 9 ) 2 - 7.870(lnE 9 ) + 0.565. (5) 

From this figure it can be concluded also that the small variation of the 
nuclear mass causes the small changes in E ff also. When the nuclear mass 
increases it causes the large changes in E ff (C; LiH; CsH) [30]. Moreover as can 
be seen from Fig. 3 in last case the empirical dependence ln<9E g /<9M ~ f(lnE s ) 
is very close to the linear one and in ordinary coordinate system it has a next 
expression: y/dE g /dM = E g /KI ° E (dE g /dM = E^ 0525 ). By the way it should 
be noted that at the large changes E g result in the changes of the force constants 
at the isotope effect from the large variation of nuclear mass. Observable in Fig. 
3 rather large scattering data in the rate of change E 9 on the isotope mass in 
the first step cause the different degree of the isotope-induced-disorder. The 
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last effect, as is well-known (see e.g. [4]), due the different magnitude of the 
scattering potential at the isotope substitution. 

We conclude the discussion of the band-gap shifts with isotope mass by ob- 
serving that the effects are exceedingly small in the most cases (see Table 2) 
and most likely will not have any consequences any semiconductors technolog- 
ical applications. However, this exercise is an excellent demonstration of the 
advanced state of our quantativc theoretical understanding of the subtle effects 
of temperature, pressure and isotope mass on the electronic band structure. 

Thanks are due to Dr. M. Segall for a critical reading of the manuscript. 
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1.0.3 Figure Captions. 

Fig. 1. Mirror reflection spectra of crystals: LiH, curve 1; LiH 2 ,D 1 _ :I ,, curve 
2; and LiD, curve 3 at 4.2 K. Light source without crystal, curve 4. Spectral 
resolution of instrument indicated in the diagram. All curves have the same 
energy scale. 

Fig. 2. Dependence of the interband transition energy E 9 in mixed LiH^Di-a, 
crystals on the hydrogen concentration x. The stright dashed line is the linear 
dependence E g = f(x) in the virtual crystal model. The solid line corresponds 
to calculations using the polynom of second degree [2] Points derived from re- 
flection spectra indicated by crosses, and those from luminescence spectra by 
triangles. 

Fig. 3. The dependence of ln(9E 9 /<9M) ~f(lnE g ); points are experimental 
data from Table 2 and continuous line - calculation on the fotmulae (5). 



Tabic 1 .Values of the energy of maxima in exciton reflection spectra of pure 
and mixed crystals at 2K, and energies of exciton binding Ef, band-to-band 
transitions E 9 . 



Energy, meV 


LiH 


LiH0.82D0.i8 


LiH0.40D0.60 


LiD 


Li B H (78K) 


Eis 


4950 


4967 


5003 


5043 


4939 


E2s 


4982 


5001 


5039 


5082 


4970 


E b 


42 


45 


48 


52 


41 


E 9 


4992 


5012 


5051 


5095 


4980 



Table 2. Values of the coefficients dE/dM (meV) and energies Of the band- 
to-band transitions E g (eV) according to indicated references. 
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13 C ^12 C 


1 4 fi n m 


k 41 25 n 01 


7 LiH^ 7 LiD 


103 \2 4 


4 992-^5 095 \2 41 


7 T,iH^ 6 T,TH 


1Z [OJ 


4 980 \1 41 




fif) 291 


4 440 291 


30 Si^ 28 Si 


2 [8] 


3.652 [8] 


68 Zn0^ 64 ZnO 


0.372 [16,20] 


3.400 [16,20] 


Zn 18 0->Zn 16 


3.533 [16,20] 


3.400 [16,20] 


65 CuCl^ 63 CuCl 


-0.076 [11,12] 


3.220 [11,12] 


Cu 37 Cl^Cu 35 Cl 


0.364 [11,12] 


3.220 [11,12] 


Cd 34 S^Cd 32 S 


0.370 [21] 


2.580 [21] 


110 CdS^ 116 CdS 


0.040^-0.068 [22] 


2.580 [22] 


Ci4 8 0->Cu£ 6 


1.116 [17-19] 


2.151 [17-19] 


71 GaAs^ 69 GaAs 


0.39 [15] 


1.53 [15] 


76 Ge ^72 Ge 


0.225 [13] 


1.53 [13] 


76 ^73^70 Ge 


0.37 [24-28] 


0.74 [24-28] 
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